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Abstract 

Considering the confinement of massless u and d quarks and massive s quarks 
in an axially symmetric quadrupole shaped deformable bag, the stability prop¬ 
erties of color-singlet strangelets at finite temperatures are studied for the 
baryon number A < 100. It is found that, in general, at each A the en¬ 
ergy of the strangelet is the lowest for the spherical shape. However, it is 
interesting to find shell structures for the deformed shapes at some new A 
values, not found for the spherical shape. This should be of relevance for the 
experimental search of strangelets. The color projected (singlet) calculation 
shows that even for the deformed shapes the shell structures vanish only for 
the temperature, T ^ 30 MeV. On the other hand, from a somewhat crude 
estimate we also find that at a fixed entropy per baryon (instead of a fixed 

temperature) of the order of unity or higher these shell structures disappear. 
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I. INTRODUCTION 


The stability properties of very small strangelets (baryon number A < 100) were stud¬ 
ied recently by Gilson and Jaffe [[TJ and others j2| in the independent particle spherical shell 
model picture where massless u and d quarks, and massless or massive s quarks are con¬ 
sidered to be confined in a spherical bag with a bag pressure constant, B. In view of the 
hot environment of ultrarelativistic heavy-ion collisions where strangelets may be formed, 
we have recently extended the work of Ref. [1J to finite temperatures using the method 
of statistical mechanics. At a finite temperature (T > 0), for a given baryon number, A, a 
partition function is constructed using the eigen-energies of the non-interacting (3A) quarks 
confined in the bag. Then the free energy of the system is minimized with respect to the 
radius of the bag leading to the minimum energy of the system at that T. For T = 0 we 
actually consider T = 1 MeV and find ||.fjj that the results are similar to those in Ref. [|lj . 
At high temperatures (T ^ 10 MeV) we construct explicitly, following projection techniques 
HQ, a color-singlet partition function. At a given T and A the color projection (CP) leads 
to lowering of the energy as compared to the unprojected case and when there are shell 
closures, the shell pockets become deeper. After CP some of the shell closure positions also 
get shifted to other A values. Using B 1 ^ 4 = 145 MeV with massless u and d quarks, and 
massive s quarks ( m s = 150 MeV) we found in Ref. |3J] that in the bulk limit, A > 100, 
CP has little effect and results are similar to those of the liquid drop model calculations of 
Jensen and Madsen [|]. It is also found in Ref. f|] that for T ^ 30 MeV shell structures 
melt away. 

Now, coming back to the main topic of shape deformation, we know that many nuclei are 
deformed in the ground state, and this happens simply due to the effect of shell structures 
(quanta! effect). Therefore, it should be worthwhile to investigate the stability properties of 
strangelets allowing for the deformation degrees of freedom at T = 0 as well as at T > 0. 
Using the same formalism as earlier [Q these calculations can be carried out if the quark 
eigen-modes in the deformed bag are known. Actually following the formalism developed by 
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Viollier, Chin and Kerman j7|J to compute the single-particle energies of quarks and gluons 
in a spheroidal bag, we have already made a study 0 of the thermodynamics of a deformed 
bag with A = 1 considering two flavors of quarks u and d, and gluons. In the present work 
the considered values of temperatures are very low (T — 1 and 20 MeV) and so like in Refs. 
|]|||] here too we have not included gluons in our calculations. 


II. FORMALISM 


The formalism is essentially that of Ref. [|J except that now the input single-particle 
energies depend on the quadrupole deformation parameter which is an additional degree of 
freedom. The computations being quite time consuming we have considered only two values 
of the temperature T = 1 and 20 MeV. 

Introducing a deformation parameter, D, the surface of a spheroidal bag, in the body- 
fixed frame, is given by 


Dx 2 + Dy 2 + irV = R 0 2 , 


( 1 ) 


where R 0 is 
is clear that 
respectively, 
spherical. A 


the radius of a sphere with the same equivalent volume. ^Frorri Eq. (jl) it 
the semi-major and semi-minor axes of the spheroid are DR 0 and RqD -1 ^ 2 , 
with a volume clearly independent of D. For D = 1 the shape is obviously 
radial coordinate at the surface of the spheroid is given by 

-1/2 


R(9) = Ro 


a(D) + J-b{D)P 2 ( cos9) 


( 2 ) 


with 


and 


a(D) = i (2 D + D~ 2 ) , 

m = -JI(d - D- 2 ) 


( 3 ) 


( 4 ) 
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where P2(cos0) is the standard Legendre polynomial of order 2. The details of the cal¬ 
culations of the eigenmodes for the quarks confined in the cavity described by Eq. (HD or 
Eq. (^) are as given in Refs. The problem is reduced to the diagonalization of a real 

non-symmetric matrix for a fixed parity. The eigenvalues are functions of Ro and D where 
R 0 enters actually as hc/R 0 with he = 197.32 MeV fm. 

The color-singlet partition function is given by |4]|| 


Z C (T,R 0 ,D) = f d^(g) Tr (u(g) e~ A / T ) , 
7su(3) v 7 


( 5 ) 


where H denotes the Hamiltonian of the system and 


with Nq = 3 for SU(3) color group. The chemical potential (/i 
function can be written as [Q 


'SU(Nc) 


9k S 2 


( 6 ) 


q = — Hq) dependent partition 


Z c (T,Ro,D,n q ) = f dn(g) e° , 
7su(3) 

with 


m 


e = \ E 1» (l + /f + Vi C08 »,) + 1E1“ (l + St + Vi cos 6 , 

+ ~ E^ 11 I* + ft + 2/^COs(^i + 62 )} + X E ' n (l + Si + 2fi COS#,) 


5E'n(l + /f + 2/jco 3 9 2 ) + iE'o 


1 + fa + ^/q 008(01 + 02 ) 


( 8 ) 


where and /| = 

I 11 Eq. (HD (3 = 1/T is the usual inverse temperature with T taken in the units of energy 
(here MeV) and e" represents the Rq and D dependent quark (antiquark) single-particle 
energies with a denoting the single-particle state labels (spherical or deformed). The upper 
limit of the summation a depends upon the value of A and T, such that further inclusion of 
states in the summation does not affect the results. In this sum the summation over quark 
flavors u, d , and s is also understood. The color unprojected (CUP) partition function is 
simply given by 
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Zu(/3, i?o, D, n q ) 


( 9 ) 
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where now 


0 = ^ [ln{l+ e -^( £ ?-^)} + In{1 + e -0( e ?+**«)}' 

Q! 

As usual the particle number for a given flavour is calculated from 


( 10 ) 


N, = r/-(ln2) , (11) 

U flq 

where 2 is given by Eq. (|7) or Eq. (9), as the case may be. The baryon number A is fixed 
by adjusting the value of the quark chemical potential such that the excess number of q over 
q is 3 A: 


(N u - N a ) + (N d - N d ) + (N s - Ns) = 3A . (12) 

The free energy of the system is given by 

F(T,Ro,D) = —T In Z + 3qi q A + BV . (13) 

With the above definition of free energy, the condition of stability is given by 

p --(w F(r ^ D> ) Tu r°' (14) 

This way at a given T and D, the bag radius parameter i?o is determined. We may call Rq 
a radius parameter in the sense that for a deformed bag it is not actually the radius. In 
practice we minimize the free energy with respect to Rq, keeping at each step the number 
Eq. (|l2|) satisfied. Then the energy of a strangelet is computed from 

E{T,Rq,D) = T 2 ^(\nZ) + 3 N A + BV , (15) 

where BV is the bag volume energy. 
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III. RESULTS AND DISCUSSIONS 


After solving Eqs. (0) and (|14|) selfconsistently for a given A, T and D the energy 
(mass) of a strangelet is computed according to Eq. As in our earlier works [|.[jj, here, 

too, we have taken B l ^ A = 145 MeV, massless u and d, and massive s quarks, (m s = 150 
MeV). The variation of the energy per baryon (E/A) as a function of A is plotted in Fig. 
la for T = 1 MeV (practically the ground state) for spherical and prolate shapes with D = 
1.0, 1.2, 1.4, 1.6, 1.8, and 2.0. A similar plot is shown in Fig. lb for D = 1.0, 0.8 and 0.6; 
for D < 1 the shape is oblate. From Eq. (]Ij) it is clear that the ratio of the semi-major axis 
(a) to semi-minor axis ( b ) of the spheroid is D 3//2 from which (just to give an idea of the 
extent of deformation) it also follows that approximately a = 0.5 b, 2b and 3 b for D = 0.6, 
1.6 and 2.0, respectively. These a and b should not be confused with those defined through 
Eqs. @ and (f|). From Figs, la and lb it is clear that, in general, the shell structure is 
persisting even if the shape of the bag (strangelet) is deformed, though energitically at every 
A the spherical shape leads to the lowest energy configuration. The shell positions at various 
baryon numbers (A s h) are shown more explicitly in Table I. One new thing we find is that 
with the increase of deformation (prolate as well as oblate) the shell structure at A = 18 is 
vanishing and, on the contrary, the shell structure at A = 14 is becoming more pronounced. 

In Figs. 2a and 2b we display the variation of the strangeness (actual value is negative) 
as a function of A for the prolate and oblate shapes, respectively at T = 1 MeV. In each 
figure the D = 1 (spherical) case is presented for the sake of comparison. In the spherical 
case we have steplike structure with well-known loading - unloading characteristics near A 
— 20 [2J. With the increase of deformation the curves are becoming smoother, though still 
some oscillation persists, as expected, at T = 0. 

In Fig. 3 we again show the plot of E/Avs A for various prolate shapes, like in Fig. la, 
but now at T = 20 MeV when the color projected (singlet) and color unprojected (CUP) 
curves become distinguishable. For the CUP case there is hardly any shell structure except 
near A = 6 for D = 1.0. As can be seen from Table II, this is actually at A = 5 with a 
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small pocket of about 1.0 MeV with the number of s-quark, n s = 2 (see strangeness vs A 
plot in Fig. 4). On the other hand, for CP case and 11=1.0-1.6, the shell structure appears 
only at A = 4 with n s = 0. For A > 4 the strangeness number becomes non-zero, so much 
so that even at A = 5 and D — 1, n s = 2.3. As Fig. 4 shows, for A > 10 or so, the S' vs A 
curve becomes almost a straight line ( S/A ~ 0.5-0.6) with the strangeness content slightly 
increasing with the increas of deformation. One striking thing in Fig. 3 for the CP case is 
the development of a shell structure at A = 15 with the increase of deformation. As can be 
seen from Table II, this shell pocket at A = 15 is 0.6 MeV with number of s-quark, n s ~ 7. 

^Frorn the discussions so far, it is clear that the color-singlet constraint leads to more 
pronounced shell structures. For T = 20 MeV and 71=1.0, 1.4 and 2.0 we demonstrate it 
in Fig. 5 through a plot of a type of color-corelation energy A E = E CVP jA — E CP /A vs 
A. It clearly shows the magnified effect of color projection near the shell closures. In this 
figure the first high peak is exactly at A = 4 whereas the second one is at A = 14 with A = 
13 and 15 lying quite nearby. It is interesting to notice that, though for D = 1.0 there is a 
clear-cut shell pocket at A = 14 in the E/A vs A plot, the effect of the gain in energy due 
to color projection is getting highly focussed in the present plot. On the whole for A > 6 
the effect of color projection seems somewhat stronger for D > 1 but again for D > 1.2 and 
A > 20 it shows a kind of saturation as far as the large values of D are concerned. 

Finally we want to discuss the effect of deformation vis-a-vis mass of the s-quark. Like 
in our previous work ||, here, too, we have considered two values, rn s = 0 and 150 MeV. 
The E/A vs A plot at T = 1 MeV is shown in Fig. 6 for m s = 0 (solid curve) and m s = 
150 MeV (dashed) for D = 1.0 (lower curve), 1.4 (middle) and 2.0 (upper). Again at each 
A the energy is the lowest for the spherical shape. The positions and the magnitudes of 
the shell structures are also enumerated in Table III for m s = 0 (see Table I for rn s = 150 
MeV). For m s = 0 a good (deep) shell structure at A = 6 persists even with the increase 
of prolate deformation (oblate not considered) unlike in the case of massive s— quark. Also 
at A =18 a deeper shell structure develops with the increase of the prolate deformation in 
contrast to the opposite trend for rn s = 150 MeV at the same baryon number, but rather 
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similar to that at A = 14 when m s = 150 MeV. 

For the m s = 0 case at T = 20 MeV only the CP case is considered and the E/A vs 
A plot for both m s = 0 (solid) and 150 MeV (dashed) is displayed in Fig. 7 for D = 1.0 
(lower), 1.4 (middle) and 2.0 (upper). For m s = 150 MeV and D = 1.0 and 1.4 the only 
shell structure seen is at A = 4 with n s = 0, whereas for D = 2.0 a shell structure has 
developed at A = 15 (see Table II). It is rather surprising to see that with m s = 0 the shell 
structure at A = 6 is 27.11 MeV, 22.34 MeV and 16.31 MeV deep for D = 1.0, 1.4 and 2.0, 
respectively which are even slightly deeper than that at T — 1 MeV. Though this feature 
appears very interesting, it can not be given much importance as m s = 0 is unphysical. But 
it does highlight that the quanta! features of the strangelets, particularly for A < 50, are 
quite sensitive to the mass of the s— quark. 

Before presenting finally our conclusions in the next section we wish to make some 
comments, in view of the work of Ref. [||], on the shell structure features if the E/A vs A 
plot is made at a fixed entropy per baryon (s). ^Frorri our color projected numbers at D=1 
and T — 1,10,20,30 MeV we have computed s = (E/A — F/A)/T at these temperatures. 
Then a crude graphical interpolation analysis shows that for s ^ 1.0 there appear no shell 
structures. It seems somewhat surprizing but a very interesting result in view of what we 
have discussed so far at fixed temperatures. 

IV. CONCLUSIONS 

We have investigated the stability properties of color-singlet strangelets at finite temper¬ 
ature within an axially symmetric quadruploe shape deformable MIT bag model. We find 
that at each baryon number A( < 100) energetically a spherical shape is preferred at zero 
as well as at a finite temperatures. However, with the increase of deformation we find that 
the shell closures at some A for the spherical shape disappear and some new ones appear at 
different A values. This should have important implications in the experimental search of 
strangelets. The maximum deformed (D = 2.0) prolate shape is considered such that the 


ratio of the semi-major to semi-minor axis is about 3:1. On the considered oblate side (D = 
0.6) this ratio is 1:2. At a finite temperature, T — 20 MeV in the CP calculation the shell 
closure occures at A = 4 for D = 1.0, 1.2 and 1.6 which is actually a lump of u and d quarks 
only. We may add that this point was missed in our earlier work of Ref. |4|]. However, at D 
= 2.0 and T — 20 MeV a shell structure is found at A = 15 with n s = 7. 

Considering massless s quarks (m s = 0), we get, in general, results similar to as discussed 
above (see Table III). However, for this case one interesting result is found that is at T = 20 
MeV the shell structure at A = 6, for D = 1.0, 1.4 and 2.0, is even slightly more pronounced 
than that at T = 1 MeV. 

It appears that in an E/A vs A plot at a fixed entropy per baryon (s) there are no shell 
structures even for low values of s ~ 1.0 and A < 30. A detailed shell model calculation 
within the present approach for s <1.0 and A < 30 should be interesting for a proper 
analysis of the shell structures. 
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TABLES 


TABLE I. Shell positions at various baryon numbers (A s h) and the value of the pocket depth 
(E s h) defined as E s h = (E/Ab — E/A s h) where Ab > A s h and E/Ab has got the highest value. 
These are shown for T = 1 MeV at a few values of the deformation parameter D. 
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TABLE II. Same as Table I for T = 20 MeV. Now color unprojected (CUP) and color-singlet 


(CP) both the cases are shown. 


D 

^4sh 

A B 

E sh (MeV) 

Case 

1.0 

5 

6 

0.9 

CUP 

1.2 

5 

6 

0.2 

CUP 

1.0 

4 

7 

6.0 

CP 

1.2 

4 

6 

4.7 

CP 

1.6 

4 

6 

0.2 

CP 

2.0 

15 

17 

0.6 

CP 
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TABLE III. Same as Table I 



-k At A = 18 a shell has developed with the ii 
= 24 and 54. 


1 


m s =0 for CP at T = 20 MeV. 


Ab 

E sh (MeV) 

9 

26.43 

19 

0.25 

28 

3.75 

50 

4.50 

59 

1.0 

100 

4.86 

8 

18.55 

20 

0.98 

26 

0.89 

51 

3.59 

100 

4.14 

8 

13.66 

21 

4.11 

51 

3.52 

100 

2.65 


of deformation and has vanished at A 

























FIG. la. Baryon energy per particle (E/A) as a fucnction of A for various values of the 
deformation parameter (.D=1.0, 1.2, 1.4, 1.6, 1.8 and 2.0) with spherical and prolate shape 
at T— 1 MeV. The results are same for both colour projected (CP) and unprojected (CUP) 
cases at T =1 MeV. 
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FIG. 4. Same as FIG. 2a with T=20 MeV for the CP (continuous lines) as well as CUP 
(dashed lines) with D = 1.0, 1.4 and 2.0 (Down to up). 







FIG. 5. The color-correlation energy defined as E CIJP / A—E CP / A at T = 20 MeV as function 
of A for D =1.0, 1.4 and 2.0. 








FIG. 6. Same as FIG. 1 for T — 1 MeV with m s = 0 (continuous lines with D = 1 (lower), 
1.4 (middle) and 2.0 (upper)) and 150 MeV (dashed lines with D =1 (lower), 1.4 (middle) 
and 2.0 (upper)) in the CP scheme. 












